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Abstract 



Most of the Ultra High Energy Cosmic Ray (UHECR) experiments and projects 
(HiRes, AUGER, TA, EUSO, TUS,...) use air fluorescence to detect and measure 
extensive air showers (EAS). The precise knowledge of the Fluorescence Light Yield 
(FLY) is of paramount importance for the reconstruction of UHECR. The MACFLY 
- Measurement of Air Cherenkov and Fluorescence Light Yield - experiment has been 
designed to perform such FLY measurements. In this paper we will present the re- 
sults of FLY in the 290-440 nm wavelength range for dry air and pure nitrogen, both 
excited by electrons with energy of 1.5 MeV, 20 GeV and 50 GeV. The experiment 
uses a 90 Sr radioactive source for low energy measurement and a CERN SPS e~ 
beam for high energy. We find that the FLY is proportional to the deposited energy 
(Ed) in the gas and we show that the air fluorescence properties remain constant 
independently of the electron energy. At the reference point: atmospheric dry air at 
1013 hPa and 23°C, the ratio FLY/Ed=17.6 photon/MeV with a systematic error 
of 13.2%. 
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1 Introduction 



The physics of the Ultra High Energy Cosmic Rays (UHECR) is a challeng- 
ing search in the field of cosmic ray physics. There is a special interest to 
measure the energy spectrum in the region of 10 20 eV where the GZK cutoff 
is predicted (lj). For more than 40 years, many experiments looked for such 
UHECR events but only a few dozen of them have been recorded so far. The 
lack of statistics and the discrepancies between the results of the two main 
experiments AGASA (@) and HIRES (@), requires new experiments to solve 
the UHECR and the GZK puzzle. In the coming years, with the deployment 
of ground based detectors: the Pierre Auger Observatory (0) and the Tele- 
scope Array (5), this number of events will increase dramatically reaching 
several hundreds. In the longer term, space-based experiments (OWL, EUSO, 
TUS) (@) are foreseen, with the goal of reaching a thousand events per year. 

For most of the past, present and future experiments, the detection technique 
relies, at least, on the measurement of the air fluorescence produced by the 
incoming cosmic ray showering in the atmosphere. This light, observed in the 
near-UV region, is induced by the de-excitation of the molecule of air (mainly 
N2) occurring along the development of the extensive air showers (EAS). This 
Fluorescence Light Yield (FLY) is rather weak (~ 4 photons/m emitted in 
47r), and depends significantly on the atmospheric conditions (temperature, 
pressure, gas components: H 2 0, Ar, ...). 

Since the 1960's, when it was first proposed to use air fluorescence for UHECR's 
detection, the FLY has been studied. In 1967, A. N. Bunner summarized the 
existing data in his thesis (0). At that time the uncertainty on the FLY was 
~30 % ! Later, the discrepancy between the AGASA and HiRes experiments 
led the UHECR community to pursue their efforts on the FLY measurement. 
In 1996, in the context of HiRes and Telescope Array, a new set of results (j8|) 
was published. More recently, M. Nagano et al. released new measurements 
on the FLY with systematic uncertainty of ~13.2 % (0). 

A better understanding of the FLY, with a goal of setting the systematic 
uncertainties below 10 %, is now needed for the UHECR experiments. Up 
to now only the pressure dependence was really measured. To improve our 
knowledge of air FLY and its behavior with respect to pressure, temperature, 
humidity, electron energy, shower age, etc, new experiments are needed. An 



overview of the current experiments can be found elsewhere (1 1 Of ) 



This paper is organized as follows: Section 2 describes the fluorescence mecha- 
nisms and presents a predictive model. In Section 3 the experimental setup is 
described. In Section 4 the data taking and signal processing are discussed and 
in Section 5 the calibration and systematic studies are presented. Finally our 
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results are described in Section 6 and a comparison with other experiments 
are discussed in Section 7. 



2 Air fluorescence mechanisms and production model 



The Fluorescence Light Yield depends on two competitive processes: excitation 
and de-excitation of the air molecules. 

In an EAS, the air is excited by the high energy charged particles of the shower. 
When such a particle traverses the air, it ionizes air and produces secondary 
low energy electrons which will excite N 2 molecules in the low energy states 
(10-20 eV), which fluoresce. The energy (and hence the wavelength) of the 
fluorescence photons corresponds to the energy difference between two excited 
states of an air molecule: the spectrum varies in a large band, from UV to IR. 
In the UHECR experiments, the photons are detected in the near UV band, 
which ranges in the 290-440 nm window in MACFLY. 

A molecular excited state can be defined by its molecular orbital and by 
its vibrational state. The air fluorescence spectrum at atmospheric pressure 
is a band spectrum dominated by IN system of N 2 and 2P system of N 2 . 
IN system corresponds to all the transitions between N 2 orbitals B 2 T^ and 
X 2 E+, and 2P system to all the transitions between N 2 orbitals C 3 U U and 



B 3 U g duJ). 



During the de-excitation, there is a contest between the radiative and the non- 
radiative processes, which are both characterized by a time scale. The first one 
may produce photons in the experimental UV band while the second process 
dissipates the energy via thermal processes. 

The mean lifetime r e of an excited state in non-isolated conditions verifies: 

1 



where the first term r Pn is the mean lifetime of the isolated excited state. For 



instance the typical values of the 2P system are ~40 ns (Il2l ) 



The non-radiative de-excitation of the molecules (quenching) comes from col- 
lisions of the excited molecules with other air molecules. The collision times 
t* c of the excited state e with the molecule of type i (i — N 2 , 2 , H 2 0, ...) are 
given by the kinematic theory of gases: 



nkTm N2 I 2m" 



T e c pi 4cr i ]] mN2+m . 
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where T is the gas temperature, P % is the partial pressure of molecules i, k 
is the Boltzmann constant, m^ 2 and rrii are the molecular mass of N 2 and 
of molecule i, and a l e is the cross section of the collision. The t\ c have been 
measured since a long time and their values vary inversely with the pressure 
P. At atmospheric pressure the quenching time of the 2P system is a few ns. 

One often expresses the relationship between the lifetime and the pressure via 
the formula: 

11 P 

7 = -0 + F ). (2) 

P' being the pressure for which the quenching processes have a collision time 
equal to the radiative de-excitations lifetime. 



In our model, the quenching processes are characterized by the quenching 
factors k l e , independent of the pressure, defined by: k\ = , 4 Jj J- 

-»/ 7Trtl 111 » 



/ mjv 9 +rrii 



2m, 



Assuming that a % e do not vary with temperature, implies that the k l e vary solely 
as the inverse square root of the temperature: k % e {T) = ^e(T )\/W- The values 
of ^e(To) are gi ven f° r a reference temperature T . 

Thus, t 1 can be re-written as: 



T ^ (jP,T) "Px fxfc«(T) • (3) 



For any excited states, e, of a gas at temperature T and pressure P containing 
a fraction / l of molecules i, the probability Q to have a radiative de-excitation 
is: 



Q 



'e 



'e 



1 

l + PTe Efk 



To 

V(To)V T 



For all the air fluorescence lines, at wavelength A, corresponding to an excited 



state e, one can express the 



FLY 



ratio of an arbitrary admixture of N2, O2 



and H 2 0, as a function of pressure P and temperature T, via the formula: 



FLY 



Xn 2 



FLY 



(P,T) = 



E d 



N 2 



X,P=0 



l + P-r eo Efk 



(4) 



where ^^|^p =0 is the limit of the FLY ratio for pure nitrogen at P=0 hPa 
(where no quenching effect by molecular collisions is expected) and xn 2 is the 
mass fraction of nitrogen in the gas admixture (i = 2 , N 2 , H 2 0). 
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In our model we assume that the fraction of the deposited energy in the gas 
converted in nitrogen excitation is Xn 2 - The parameters entering equation 4: 
^|^ 2 P=0 , r eo , kgfa, k°fa and are given for 24 wavelengths (A) cor- 

responding to 5 excited states (e) at the reference temperature T = 20°C 
(293.15 K). The values of these parameters come from both our MACFLY 
measurements and from already published results: (0) (0) (0) for the spectra 



(A) and (jl2j) (Il3l ) for the lifetimes (r eo ) and the quenching factors (fc 



In this paper, we will show experimental results for pure nitrogen (iV 2 ) and for 
our experimental dry air (DA) — an admixture 80%(A^ 2 )-20%(O2) — as a func- 
tion of the pressure. From this it will be then possible to derive the behavior 
of the atmospheric dry air (ADA) — an admixture 78.08%(iV r 2)-0.93%( J 4r)- 
20.99%(O2) — for which we assume the FLY for argon to be identical to the 
FLY of nitrogen due to its catalysis effect on nitrogen fluorescence (flil ). Then 
the ADA becomes in our model — an admixture 79%(A^)-21%(02). 

Using this, we formulate the general FLY of any humid "atmospheric air" 
(AA), i.e. any admixture composed of (1 — /j) dry air (ADA) and \i water 
vapor, where \i is the molecular fraction of water in the admixture (/x = f H2 ° 
as defined in equation 3). We get: 



FLY 



E r] 



\ FLY 
{P,T,fi) = —- 
a h, d 



ADA 1 _ 18 



x j - 29 ^ N , (5) 

y<T 0) p e(To) t 



\,P=0 



where 



pi AD A t j piH 2 1 (a) 

e " r eo ■ (0.79^ + 0.21^) Q e r eo ■ k?'° ' 1 ' 



We assume that the energy fraction transferred to ADA molecules is propor- 
tional to the mass fraction of ADA in the admixture rather than the molecular 
fraction (1— fx). The numerator (1— ||/x) represents this mass fraction. It is a 

simplification of the exact formula (i_^29+jl-i8 ^ or A* as ^ * s ^ ne case m 
the atmospheric air. 



3 Experimental setup 



The MACFLYuD experiment is twofold, it has been designed to measure 
both the FLY induced by single electron track and the FLY produced by a 
high energy electromagnetic shower developing in the air. It is composed of 

2 Measurement of Air Cherenkov and Fluorescence Light Yield 
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two devices: MF1 which is used for single track FLY measurements and MF2 
which measures the fluorescence produced by an electromagnetic shower. In 
this paper we will present only the results obtained with MF1 using both a (3~ 
radioactive source (Strontium 90) and an electron test beam (CERN/SPS-X5 
line). The results of the shower FLY with MF2 will be presented elsewhere (fl6l ). 

The MF1 experiment is composed of a pressurized chamber (with pressure 
from to 1200 hPa) containing the gas, equipped with a trigger system (see 
Figure 1). This device has an internal cylindrical chamber, 150 mm in diameter 
and 288 mm long. An ionizing particle (E > 1 MeV) reaches the gas volume 
through the entrance window (0.25 mm of black Delrin), crosses the gas volume 
and leaves the chamber at the exit window (0.8 mm of aluminium) where it 
reaches the trigger plane. 



PMT EMI 9820QA 
Quartz window 




25 cm 



28.8 cm 



Fig. 1. Schematic view of the MF1 chamber (left) and a cut view of its optical 
system (right) showing the multilayered mirror, the quartz lens, the Winston cones, 
the BG3 filters and the PMTs. 



An optical system collects the induced fluorescence light, and focuses it on 
two UV sensitive phototubes (PMT) EMI9820QA. Each PMT receives pho- 
tons via an optical system which is composed of a multilayered mirror covering 
the chamber, a quartz lens focusing the light, a Winston cone and wide band 
filters which allow the selection of the appropriate wavelength. The multilay- 
ered mirro E was designed to have its best reflectance in the appropriate air 
fluorescence band (see Figure 2). The filters are Schott BG3 filters with a large 
band of transmittance: 290-440 nm. We have measured the FLY for two gases: 

3 From internal to external layer: 90 nm (Al) + 43 nm (S1O2) + 43 nm (Hf02) + 
43 nm (Si0 2 ) + 43 nm (Hf0 2 ) 
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Fig. 2. MF1 multilayered mirror reflectivity coefficient as a function of wavelength 
(upper curve) and BG3 UV filter transmission coefficient (lower curve). 




Scintillator 



rigger PMT 



Electrons 



-PMT 



Fig. 3. MF1 setup in its laboratory configuration. The radioactive source is embed- 
ded in a shielded box, and the trigger is mounted on the device itself. 



pure nitrogen and dry air (A^-80%, 02-20%). The MF1 experiment was tested 
in two different setup configurations corresponding to low and high energy. For 
the low energy runs (labeled MFl-lab) a f3~ source ( 90 SV) was used and the 
trigger system was integrated in the overall system (see Figure 3). For the 
high energy runs (labeled MFl-beam) the CERN SPS XP5 beam facility was 
used. In this later case the source and the trigger system are removed from the 
main tank, and external systems are installed along the beam line: a position 
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sensitive (er ~ 0.5 mm) X-Y delay chamber upstream from the MF1 device 
to record both the horizontal and the vertical beam position and two double 
plane trigger systems before and after the MF1 setup. The electron beam was 
a pulsed beam of about 10 000 electrons per spill (4.8 s duration) every 16.8 s 
with a beam spot of about 4x7 mm 2 . The electron energy was selected in the 
range between 10 and 100 GeV. 



4 Data taking and signal fitting 



The data was collected on an event per event basis. For every event the PMT 
signals (MF1 and trigger counters) were recorded by a QADC (CAEN-V792) 
which integrated the charge during a gate of 100 ns. The beam position was 
also recorded for the MFl-beam configuration. During a run, which consists 
of typically 10 6 triggers, we recorded two kinds of events: the beam events 
(BE), used for the FLY measurement, and the random events (RE), when 
no beam is present, for the background estimations and studies. In MFl-lab 
configuration BE are triggered when an electron reaches the trigger plane after 
the exit window. In MFl-beam configuration, the electron should reach two 
trigger planes before and after the MF1 chamber. In the analysis, we also 
suppressed BE not detected in the central region of the delay chamber. 

The FLY is rather weak (~4 ph/m) and most of the produced photons are 
lost in the chamber. Among those reaching the PMT photocathode, some 
are converted in photoelectron (pe). Overall, the mean number of detected 
photoelectron is about 0.01 pe/event. To extract the mean Detected Light 
(DL), we performed a fit of the data by a function describing the expected 
signal; for every event, a PMT detects an integer number of photon: 0, 1, 2 
or more converting them into photoelectrons. So the overall PMT spectrum 
is described by a weighted sum of individual photoelectron contributions, the 
weights following a Poisson distribution. 

The zero photoelectron contribution, the pedestal, is measured experimentally. 
The single photoelectron contribution is described by two functions: one de- 
scribing the standard photoelectron multiplication at the first dynode, the 
other corresponds to the photoelectron being inelastically back-scattered on 
the first dynode (17). The first function is described by a WeibulQ distribu- 



tion (Il5l ) while the second uses an exponential law. 

The Multi-photoelectron contributions (n>2) are described by gaussian distri- 
butions. The mean value \L n and the variance a 2 n of each distribution are deter- 
mined by the single photoelectron fit function parameters: /i n = n/ii + /io and 
a 2 = na\ where \i\ and o\ are the single photoelectron distribution (Weibull) 



Weibull(x) = c/b ■ (x/b) c ■ exp(—(x/b) 
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mean value and variancq_i| and fio the pedestal mean value. 




100 120 140 160 180 200 220 240 

ADC bin 



Fig. 4. PMT spectrum fit function (solid line) is a sum of several contributions: 
pedestal (black solid line), single photoelectron (solid line), multi-photoelectron 
(dotted line) and back-scattered single photoelectron (dotted dash line). 

One can see in Figure 4 the PMT raw signal as well as the different contri- 
butions entering the fit function. It is worth noting the agreement between 
the global fitted function and the data extends on more than five orders of 
magnitude. 

In the context of MACFLY data (~ 5 x 10 5 BE/run), the low level of PMT 
dark noise enables the detection of a signal as low as: 0.0002 pe/event, which 
correspond to a mean sensitivity of the MF1 apparatus of about 0.04 ph/event 
for the florescence light. Howerver other sources of PMT photoelectrons signal 
degrade this minimal sensitivity to ~0.2 ph/event. 

The detected light has several origins: fluorescence light (FDL), Cherenkov 
light (CDL) and background (Bgd). 

DL = FDL + CDL + Bgd . (7) 

The Bgd contribution is estimated from both the data from the RE triggers 
and from BE with vacuum in the chamber (no signal expected). 
The Cherenkov contribution is estimated, based on a Geant4 Monte Carlo 
simulation program, describing in details both the apparatus and the inter- 
action processes (I la ). For MFl-lab measurements, electrons do not produce 
any Cherenkov radiation. At high energy (E > 10 GeV), the Cherenkov light 
yield is important but its contribution to the raw PMT signal remains small; 

* Ml = 6.r(i + i) a nd* 1 = by/r(i + i)-(r(i + ±)) 2 
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the Cherenkov light is emitted in the forward direction where a black "light 
catcher" is installed, suppressing dramatically the reflected Cherenkov com- 
ponent. 

Figure 5 shows the DL measured and the estimation of CDL and Bgd for 
test beam measurements at CERN as a function of pressure in Dry Air and 
nitrogen. The FDL is determined by substraction. We can see that the main 
part of the DL comes from the fluorescence. The Fluorescence Light Yield in 




Pressure (hPa) 



10" 

Pressure (hPa) 



Fig. 5. Measured light in dry air (left) and in pure nitrogen (right) in milli-photo- 
electron per event (mpe/evt) as a function of the gas pressure. Triangles are raw 
measured light data (DL); dotted line is the Bgd estimation from vacuum measure- 
ments; dotted dash line is CDL simulation; stars are the FDL data (after substrac- 
tion of Bgd and CDL); dashed line is a simple FLY model fit and the solid line is 
its projection on DL. 

MF1, is determined from the measurements from both overall and background 
signals and from the estimated contribution for the Cherenkov yield: 



FLY 



DL - CDL - Bgd 

£mfi 



where the MF1 efficiency Smfi is the product of geometrical acceptance of the 
detector and of the PMT quantum efficiency {eMFi-beam = 0.00556 ph/pe). 



5 Calibration and systematics 



The geometrical acceptance is estimated by a Geant4 (Il8l ) simulation program 
which describes in detail the apparatus. It is able to perform the tracking of 
the optical photons from the production source along the electron track to 
the PMT photocathodes (Filters transmittance and Mirror reflectivity are 
included). The simulation result shows that 2.75% of the fluorescence light 
emitted isotropically in the chamber reaches one PMT photocathode. 
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Uncertainty sources 


MFl-lab 


MFl-beam 


Electron position 


5% 


6% 


Mirror 


4% 


4% 


Winston cone 


2% 


2% 


Row aluminum 


3% 


3% 


Cherenkov catcher 


0% 


1% 


Others (lens, filter ...) 


< 2 % 


< 2% 


Geometrical Acceptance 


7.5% 


8.2% 



Table 1 

Geometrical acceptance systematic uncertainties for both experimental configura- 
tions : MFl-lab (with 90 Sr source) and MFl-beam (CERN test beam) 

This geometrical acceptance varies with the setup configuration (mirror, sur- 
faces, filters,...). The optical properties of the inner surfaces are well known: 
the mirror reflectance, as function of wavelength, was measured with less than 
2% of error. Figure 2 shows the result of this measurement. The other inner 
surfaces were studied in laboratory with a calibration device at 370 nm. 

To check our Geant4 based simulation program, we compared the FLY mea- 
surements for different experimental (optical system) configurations with the 
simulation expectations. For instance, having black covered some inner sur- 
faces, the geometrical acceptance was reduced by 24% which was really close 
to the predicted simulation value of 26%. 

Table 1 summarizes the systematic relative errors on geometrical acceptance 
for both experimental configurations (MFl-lab and MFl-beam). The main 
contributions come from the electron track position uncertainty (delay cham- 
ber/MFl alignment and multiple-scattering in the chamber), from the internal 
surfaces reflective properties and from the mirror inhomogeneity. The overall 
geometrical acceptance uncertainty is estimated to 7.5% for MFl-lab and 8.2% 
for MFl-beam. 

The others sources of systematic uncertainty come from our reconstruction 
method used to extract the fluorescence signal from the data and from the 
PMT calibration. We identified three sources of errors from our reconstruc- 
tion method: PMT signal fitting procedure (DL reconstruction error), back- 
ground and Cherenkov estimations. These uncertainties vary with the setup 
configuration and depend on the signal intensity. 

Table 2 gives our estimation of the systematics uncertainty for dry air measure- 
ments with the MF1 setup for both laboratory and test beam configurations. 
The total absolute uncertainty is 13.2% for MFl-lab and 13.7% for MFl-beam. 
They are dominated by PMT calibration uncertainties. 
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Errors sources 


MFl-lab 


MFl-beam 


Geometrical Acceptance 


7.5% 


8.2% 


PMT calibration 


10% 


10% 


DL reconstruction 


4% 


3.5% 


CDL Simulation 


0% 


2% 


Bgd Measurement 


1% 


2% 


Systematic Error 


13.2% 


13.7% 



Table 2 

Systematic uncertainties on the FLY measurement for both experimental configu- 
rations: MFl-lab and MFl-beam. 

6 Results 



In our air FLY model presented in section 2, the light yield is proportional to 
the deposited energy in air (Ed). Therefore the ratio becomes independent 
of the energy lost by the incoming electron. It will only vary as a function of 
the atmospheric conditions (pressure, temperature and humidity). 

In the data used for this paper, the pressure is the only parameter under 
consideration. Figure 6 shows the measured by MF1, in the wavelength 
range 290-440 nm, as a function of the pressure, both in our experimental dry 
air (DA) and in pure nitrogen for different energies of the incoming electron: 
1.5 MeV, 20 and 50 GeV. The dotted lines correspond to our FLY model 
described in section 2 (equation 4), fitted on the MACFLY data, for which 
i = N 2 ,0 2 . We have fitted the absolute value and the pressure dependence 
of the FLY (from 3 hPa to 1100 hPa) both in our dry air (DA) and in pure 
nitrogen (in same time). The fit had 3 free parameters: the 290-440 nm integral 
value of pure nitrogen FLY and the mean value of the k^ 2 and of the 2 . 

The zero pressure extrapolation FLY obtained, respectively, for pure nitrogen 
and dry air are: 



FLY 



n 2 ply 
1959 ± 412 ph/MeV and 



p=o Ed 



DA 

= 1523 ±329 ph/MeV 

p=o 



The first integral value is used to compute the z§^}aj > =o va hies for the 24 A, 
according to the spectral published values (0) (jH) (fill ). In the same way, the 
two mean k\ values extracted from our pressure dependence measurements 
are used to compute the 10 k l e for the five exited states according the relative 



values from (Il2l ) (jl3l ). Table 3 summarizes the values of all the parameters 



used in our model. The P' are determined by the equation 6 with the fitted 



k\ and already published life times (1l2l ) (1l3l ) 



At atmospheric pressure and T = 23°C in dry atmospheric air, our FLY model 
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FLY/E d (P) in Dry Air | 
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Fig. 6. Ratio FLY ' j Ed - "FLY on deposited Energy" as a function of the pressure, 
in our experimental dry air — DA — (left) and in pure nitrogen — ./V2 — (right), 
measured for different electron energies: 1.5 MeV (W), 20 GeV (k) and 50 GeV(»). 
The dotted lines show the results of the fit of the data by our FLY model. 



fitted on MACFLY data gives: 



FLY 



ADA 



17.6 ± 2.3 ph/MeV 



(9) 



The ratio of the FLY/ for pure nitrogen compared with the one for dry air, 
is shown in figure 7 as a function of the pressure. 



X 



-!"• 



1" 



10 



10 2 



10° 

Pressure (hPa) 



Fig. 7. Ratio of the FLY/Ed m -^2 and in dry air as function of pressure. The 
dotted line is the fluorescence model fit on MF1 data 

At zero pressure this ratio tends to the ratio of nitrogen mass fraction in each 
gas (1/Xn2 = 1-286), while at atmospheric pressure the quenching effect of 
the oxygen reduces by about a factor 5.5 the FLY in pure nitrogen. 

During the data taking we measured the FLY from a "real" humid atmospheric 
air (RAA). We filled MF1 chamber with the ambient laboratory air having the 
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Excited 
state 




MACFLY Model Parameters 


Atmospheric air 


A 

(nm) 


rLl jv 3 

E d x ' p ~" 

(ph/MeV) 


(ns) 


(mslPa 1 ) 


e(20°C) 

(ms-'.Pa 1 ) 


k Hl ° 
(ms'.Pa- 1 ) 


rLl \ADA 

(ph/MeV) 


D 'ADA 

^ e[20°C] 

(hPa) 


p ,H 2 

r s [20"C ) 

(hPa) 




337.1 


227 










174 






C 11„ 


357.7 


182 










139 






380.5 


64 


37 


1.65 


43.6 




49 


25.8 


2.94 




405.9 


26 










20 








434.4 


7 










5 








315.9 


149 










114 






c 3 n 

u 


333.9 


7 










5 






353.7 


48 


40 


4.18 


54.0 


95 


37 


17.1 


2.63 


u=l 


375.6 


63 






49 


399.8 


39 










30 








427.0 


10 










8 








297.7 


31 










24 






c 3 n„ 


313.6 


52 










40 






350.0 


7 


40 


7.35 


77.0 


98 


6 


11.4 


2.55 




371.1 


23 




18 


394.3 


16 










13 








420.1 


11 










8 








296.2 


20 










15 






Cn. 


328.5 


34 










26 






346.9 


9 


40 


9.06 


113 


111 


7 


8.8 


2.25 


u=3 


367.2 


3 










6 






389.4 


5 










4 








391.4 


662 


62 


38.7 


96.8 


211 


508 


3.17 


0.76 




427.8 


259 




198 



Table 3 



Parameters used in our MACFLY model describing both a general A^C^-f/^O 
admixture as defined in equation 4 and its simplification for an atmospheric air 
(equation 5). The 24 wavelengths (grouped per excited states) and the excited states 
enter the first two columns. The other columns contain the values of the parameters. 
Among them the bolded ones are obtained by htting the MACFLY data. 



Energy 


FLY/E d (ph/MeV) 


FLY /I (ph/m) 


1.5 MeV 


17.0±2.3 


3.14±0.41 


20 GeV 


17.4±2.5 


4.22±0.61 


50 GeV 


18.2±2.5 


4.44±0.61 


All energies 


17.6±2.3 





Table 4 

Values of the FLY (FLY/E^ and FLY /I) for the 3 energy samples, at the reference 
point (P=1013 hPa, T = 23° C, \i=0) for atmospheric dry air (ADA). The last line 
shows the value of FLY/Ed when all the data are fitted together. 



following parameters: P=950 hPa, T = 23°C and relative humidity 35% (in- 
formations provided by external probes). This corresponds to an atmospheric 
air (AA) with a water vapor fraction ll = 1.05%. We compared this measure- 
ment with a result from our Dry Air at same pressure and temperature. We 
found: 

FLY DA , . 

1.11 ±0.07, (10) 



FLY RAA 

while the expected value from our model is 1.13. 
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7 Comparison with other experiments 



To compare our results with the other experiments, we have chosen a reference 
point: atmospheric dry air at P=1013 hPa (1 atm.) and T = 23°C (296 K), 
and we express the FLY per track length (FLY /I) in the usual units: photon 
per track length (ph/m). Table 4 gives values we obtained for three electron 
energies: 1.5 MeV, 20 GeV and 50 GeV. These values are corrected, according 
to our fluorescence model, to correspond to this reference point. 



In Figure 8 we show the FLY jl measurements of the corrected MACFLY data 
together with the values for others experiments as a function of the incoming 
kinetic energy. On the same figure we show the expected energy depositions in 
the gas. The experimental data come from the following experiments: former 
results of Nagano et al. (for 0.85 MeV electrons), the results of Kakimoto et 
al., at 1.4 MeV, 300 MeV, 650 MeV and 1000 MeV, and the recent results of 
the FLASH experiment (11 91 ) for 28.5 GeV electrons. 



[ FLY of dry air at 1013hPa and 23°C | 




dE/dX lost 

dE/dX deposit 

Kakimoto 

Nagano 

FLASH 

MACFLY 



0.4 



> 



o.3 e 

x 

T3 



0.2 



0.1 



10 10 10 

Kinetic Energy (MeV) 

Fig. 8. FLY per electron track length FLY /I (in photon /meter) as a function of 
electron energy measured by several experiments (see text). We also show the energy 
depositions (right axis): energy lost dE/dX (dashed curve) and the energy deposited 
dEd/dX in MF1 chamber (solid curve). 

The variation of FLY jl (left axis of figure 8) as a function of the electron 
energy is compared to the dE/dX of the electron (right axis). The beam elec- 
trons lose their energy by ionization and by producing high energy delta and 
gamma rays. The dotted line shows the total dE/dX lost by the an electron, 
calculated u sing Berger- Seltzer formula (J20|) which is used by Geant4 simula- 
tion toolkit (lla ). The solid line indicates the energy deposited (dEd/dX) in the 
fiducial gas volume of the MF1 chamber, computed with our MF1 simulation 
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program (Il5l ). The difference between the two curves dE/dX and dE^/dX 
reflects the energy carried away by the high energy 5-rays and 7-rays beyond 
the MF1 chamber. The ratio between the two scales corresponds to our FLY 
model at this reference point: FLY/ E^ = 17.6 ph/MeV. 

On a wide energy range (from MeV region to 50 GeV region) the MACFLY 
data are well described by the deposited energy distribution. At low energy our 
measurement is in agreement with results from Nagano et al. and Kakimoto 
et al., following the Berger- Seltzer curve. At high energy they also agree well 
with the FLASH result, showing the deposited energy behaviour of the FLY. 



8 Conclusions 



We have performed measurements of the dry air and pure nitrogen Fluores- 
cence Light Yields induced by single electrons of low (1.5 MeV) and high (20 
and 50 GeV) energy as a function of the gas pressure. We show that, within the 
experimental uncertainties, the FLY is proportional to the deposited energy 
in the gas, independently of the incoming electron energy. At the reference 
point: P=1013 hPa and T = 20°C, the FLY/E d = 17.6 ± 2.3 ph/MeV. 

Based on our measurements and using already published data we have pro- 
posed a model describing, in the 390-440 nm range, the fluorescence light yield 
for any air composition, as a function of the pressure, the temperature and 
the water contamination. 
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and in pure nitrogen — N 2 — (right), measured for different 
electron energies: 1.5 MeV (T), 20 GeV (A) and 50 GeV(«). 
The dotted lines show the results of the fit of the data by our 
FLY model. 

7 Ratio of the FLY/ Ed in N 2 and in dry air as function of 
pressure. The dotted line is the fluorescence model fit on MF1 
data 

8 FLY per electron track length FLY jl (in photon/meter) as a 
function of electron energy measured by several experiments 
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